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Abstract 

, In this work the Melnikov method for perturbed Hamiltonian wave equations 

is considered in order to determine possible chaotic behaviour in the systems. The 
backbone of the analysis is the multi-symplectic formulation of the unperturbed 
PDE and its further reduction to travelling waves. In the multi-symplectic approach 
two separate symplectic operators are introduced for the spatial and temporal vari- 
ables, which allow one to generalise the usual symplectic structure. The systems 
under consideration include perturbations of generalised KdV equation, nonlinear 
wave equation, Boussinesq equation. These equations are equivariant with respect 
Q | to abelian subgroups of Euclidean group. It is assumed that the external perturba- 

■ tion preserves this symmetry. Travelling wave reduction for the above-mentioned 
systems results in a four-dimensional system of ODEs, which is considered for Mel- 
nikov type chaos. As a preliminary for the calculation of a Melnikov function, 
we prove the persistence of a fixed point for the perturbed Poincare map by using 
Lyapunov- Schmidt reduction. The framework sketched will be applied to the analy- 
sis of possible chaotic behaviour of travelling wave solutions for the above-mentioned 

■ PDEs within the multi-symplectic approach. 

(N 
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^ ■ 1 Introduction 

Recently it was shown how many nonlinear PDEs can be formulated in a multi-symplectic 
form H H, [|. This formulation assigns distinct symplectic structures to the spatial and the 
temporal coordinates, thereby generalising the usual Hamiltonian formulation. By means 
of the multi-symplectic approach questions like stability of solitary waves, existence of 
generalised basic state at infinity, equivariant properties of the solutions etc. can be 
considered in a more general setting, yielding new results on these issues. 

In order to analyse chaotic behaviour of travelling wave solutions to Hamiltonian 
PDEs, Melnikov's method can be used. The problem with its direct application is due 
to the symmetry in a multi-symplectic formulation of these PDEs, which results in the 
presence of unit eigenvalue among the spectrum of the Poincare map. This complication 
is solved by means of Lyapunov-Schmidt reduction. 

Finally, we illustrate the application of the Melnikov method to the study of chaotic 
behaviour in a perturbed Korteweg-de Vries equation. 
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2 General setting 

We start to consider a multi-symplectic PDE of the form ||: 



MZ t + KZ X = VS{Z) + eS 1 {Z, x-ct), Z 



( U \ 

V 

w 
V $ J 



x G 



where M and K are constant skew-symmetric matrices on R 4 and S : R 4 — > R is suf- 
ficiently smooth ( at least twice continuously differentiable). The perturbation S± is 
assumed to be a periodic function of its last argument: Si(-,x) = Si(-,x + T) and also 
C r , r>2. 

We suppose that the system ([!]) is equivariant with respect to a one-parameter Lie 
group, whose algebra is spanned by the generator £. For the unperturbed case (e = 0), 
multi-symplectic Noether theory provides the existence of the two functionals P(Z) and 
Q(Z) such that § 

M£(Z) = VP{Z), K£(Z) = VQ(Z), (2) 
The state at infinity should satisfy 0, H|: 



VS(Z ) = aVP(Z ) + bVQ(Z 



(3) 



with P(Zq) = V and Q(Z ) = Q specified real parameters, a, b G R. 

A shape of an unperturbed solitary wave travelling at speed c, Z(x,t) = Z{x — ct), 
which is biasymptotic to this state should satisfy the equation 



Z X = J- 1 VH (Z), 



(4) 



where H (Z) = S{Z) - aP(Z) - bQ(Z), and J c = K — cM |, §. 

To study the existence of travelling waves and their chaotic behaviour, we consider 
the dynamical system ( similar consideration can be found in 0) 



—Z = f (Z) + efi{Z 1 x), Z GR 4 , 0<e<l, < x < oo, 
ax 

where f (Z) = Jj x VH (Z), and f x (Z,x) = J^S^Z^x). 
The following hypotheses are imposed on the system: 

(HI) a) f : R 4 -> R 4 is C r (r > 2); 

b) /i:K 4 x5'^ R 4 is C r (r > 2). 

The system (||) can be rewritten as a following suspended system: 

r dZ 



(5) 



dx 

d6 
< dx 



f (Z) + eh(Z, 



(6) 



with the frequency u = 2n/T. Its flow $ t e : R 4 x S 1 -»■ R 4 x S 1 is defined for all t G 
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(H2) a) The unperturbed system 



^-Z = f (Z) = J-'VHoiZ) (7) 



is Hamiltonian with energy H : M 4 — > K. 

So, we have the corresponding symplectic form 

Q(Z 1 ,Z 2 ) = (3 C Z 1 ,Z 2 ). (8) 

b) The system fl7|) zs equivariant with respect to a one-parameter symmetry group Q 

spanned by the generator q. This group Q is assumed to be either compact or a 
subgroup of affine translations. We also suppose that the perturbation preserves 
this symmetry. 

c) The system (|7|) has the family of fixed points 4>po, where po = and (j) G Q , and 

corresponding (heteroclinic) orbits Z (x) such that 

^Z (x) = f (Z (x)), (9) 
ax 

and lim^.oo Z (x) = 4>ip as well as lim^+oo Z (x) = (p 2 Po, for some 4>i,(fi 2 £ 
Q. 

(H3): fi(Z,x) = A\Z + f[x) + g(Z,x), where A\ is a linear operator, fix) = f(x + T), 
g(Z,x) is time-periodic with period T and also satisfies g{0,x) = 0, Dg(0,x) = 0. 

(H4) a) For e = the spectrum a[exp(TA)] = {1, 1, e ±XT }, X > 0, where A = Z~ 1 D 2 Hq{p g ), 
b) i. (Hamiltonian case) For e > a[exp[T(A + eAi)]] = {1, 1, e TK } , 

ii. (Dissipative case) For e > a[exp[T(A + eAi)]] = {l,A d ,e TAf }, where 
de ^ dist(\ d , \z\ = 1) < C 2 e, Ci >0,C 2 > 0. 

Next, one can define the Poincare map P e : M 4 — > R 4 as 

P e (Z) = 7n^(Z,0), (10) 

where 7Ti : IR 4 x S* 1 — > IR 4 denotes the projection onto the first factor. Equivalently, one 
can define 

P£ (Z)=ir 1 & T (Z,x ). (11) 
We rewrite the fixed point equation P e (p e ) = p e in the form: 

V e (p e ) = 0, (12) 

where V e (Z) = P e (Z) — Z, and the operator L = DV°(0) is introduced. 
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3 Main results 

Lemma l.Let (H1)-(H4) hold. For e small, there exists a unique group orbit 4>p e of 
fixed points of the perturbed Poincare map near the group orbit <ppo such that 

min - 2 Po} = 0(e). 

91, 4>2&G 

Equivalently, there is a family of periodic orbits 07 e (^) = (4>p e ,0Jx) of the perturbed system 
(dj near 07o(a;) for G Q . 

Lemma 2. For e > sufficiently small, we have a [DP e (p e )] = {1,1, e e } for the case 
(H4b i) and a [DP e (p t )] = {1, X d , e TXt } for the case (H4b ii) respectively. 

Conjecture 3. Corresponding to the eigenvalues e TA *,l there exist invariant mani- 
folds: W ss ('j e (x)) ( the strong stable manifold), W u (^ e (x)) ( the unstable manifold), and 
W /C (7 e (x)) ( the centre manifold) of p e for the Poincare map P e (Z) such that 

i. W^j^x)) and Wf^.( n f e (x)) are tangent to the eigenspaces of e TAe respectively at 
7 e , while Wf oc ( , y e (x)) is at the same point tangent to the eigenspace corresponding to unity 
eigenvalue ( double in the case of Hamiltonian perturbations) . Their global analogues are 
obtained in the usual way: 

W ss (j £ (x)) = U,< W» ( 7e (a0), n ,x 

W«Mx)) = U.>o^W2 e (7«(a?))- 1 J 

ii. they are invariant under P e (-); 

Hi. W ss (%(x)) and W w (%(x)) are C r O(e)-close to W s ( / y (x)) and W u (^q(x)) respec- 
tively. 

Lemma 4. Let (Zp u (x,xo),ujx) be orbits lying in W s,u ( , y e (x)) and originating in an 0(e)- 
neighbourhood of (Zq(-xo), 0). Then the following expressions hold with uniform validity 
in the indicated time intervals: 



Zf(x, x Q ) = Z (x - xq) + ey*(x, x ) + 0(e 2 ), x e [x , oo), 
Z?(x,x ) = Z (x - x ) + ey^(x,x ) + 0(e 2 ), x 6 (-oo,x ), 

where y s e ,u satisfy the first variational equation: 



(14) 



^ = J c 1 D 2 H (Z (x - x ))y + efi(Z (x - x Q ),ux). (15) 
We introduce Melnikov function as: 

/oo /»oo 
DH (Z (x)) ■ fi(Z (x),x + x )dx = / Q(f (Z (x)),fi(Z (x),x + x )dx. 
-oo J — oo 

(16) 

Theorem 1. Suppose M(xq) has simple zeros. Then for e > sufficiently small W /S (7 e ) 
and W u {pf e ) intersect transversely. 

This result implies via the Smale-Birkhoff theorem || the appearance of a horseshoe 
near the saddle point of the perturbed Poincare map, what results in a chaotic dynamics 
in the corresponding region of the phase space. 
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4 Example 



We consider the perturbed generalised Korteweg-de Vries equation |], ||: 

U t + AU X + aU p U x + U xxx = ef x (U, U x , U t , x), (17) 

where the perturbation is assumed to be periodic of period T in its last argument: f(x) = 
f(x + T). It can be rewritten in a multi-symplectic form as 

/ U \ 



MZ t + KZ X = VS, Z 



V 

w 
V $ J 



x G 



with 



M 



/ -1 \ 





\ 1 / 



K 



and 



S = -V 2 - -UW + -AU 2 + 
2 2 2 



/0 -1 





o \ 


1 














-1 


\ 


1 


o / 


a 




-U p+2 


(p + l)(p- 


f2) 



(19) 



(20) 



The unperturbed solution is defined as: 

U (x) = 2b + \K 2 sech 2 (^) , 
Vb(x) = -fi^sinh (^) sech 3 (^) , 
W (x) = 2a + 46(A + a6) + ^cK 2 sech 2 (^) , 
$ (x) = ^[l + tanh(f9] Q , 

Here K = \J c — A — 2a6. Perturbation can be represented in a multi-symplectic setting 



(21) 



as: 



efi(Z,x) 



( 



\ 



\ 

ef(Z,x) 








(22) 



/ 



Here f(Z,x) = f(U,U x ,—clI x ,x) with U and its derivatives being components of Z. 
Therefore Melnikov function ([16]) for the system (O) yields 



M(x ) 



V (x)f(Z (x),x + x )dx. 



(23) 



For the one-harmonic dissipative driving force this Melnikov function will have simple 
zeros 0, and therefore ona can conclude the chaotic dynamics of U. 



5 Conclusions 

We have considered a modification of Melnikov's method, which can be used for the 
analysis of chaotic behaviour of travelling wave solutions to multi-symplectic PDEs. The 
results are illustrated with the example of the perturbed KdV equation. 
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